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Introduction 

We  have  recently  [1-3]  analyzed  thermodynamic  data  for  the  solvation  of  monoatomic 
monovalent  ions  in  water  and  a  number  of  other  polar  solvents  on  the  basis  of  the  mean  spherical 
approximation  (MSA)  [4-6].  The  attractive  feature  of  the  MSA  is  that  it  provides  simple  analytical 
results  which  can  easily  be  used  by  those  analyzing  experimental  data.  An  important  departure  of 
our  work  from  that  carried  out  earlier  is  that  the  polarization  parameter  A,  which  characterizes  ion- 
solvent  interactions  is  regarded  as  an  adjustable  parameter  whose  value  is  estimated  from  the 
experimental  Gibbs  solvation  energy.  This  leads  to  the  conclusion  that  X  depends  not  only  on  the 
nature  of  the  solvent  but  also  on  the  sign  of  the  charge  on  the  ion.  Although  this  result  is  perhaps 
not  surprising,  it  emphasizes  that  a  model  which  considers  only  point  charge-point  dipole 
interactions  cannot  adequately  describe  experimental  results.  In  other  words,  the  statistical 
mechanical  description  of  an  electrolyte  solution  must  include  consideration  of  higher  order 
moments,  and  the  polarizability  of  the  solvent.  These  features  of  the  system  are  incorporated  in  a 
'sticky'  parameter  [7]  which  reflects  the  fact  that  anions  and  cations  bind  differently  to  the  solvent 
Anions  dismpt  water  structure  very  little,  for  example,  whereas  cations  which  are  attracted  to  the 
electronegative  oxygen  atom  cause  local  disruption  of  the  solvent.  It  can  be  shown  [7]  that  the 
effect  of  this  sticky  interaction  is  to  change  the  value  of  the  parameter,  X,  so  that  instead  of  a  single 
value  of  X,  which  is  derived  for  a  mixture  of  ions  in  a  solvent  of  hard  dipoles,  we  have  different 
values  for  each  ionic  species,  Xj.  All  deviations  from  the  MSA  will  be  lumped  into  this  effective 
polarization  parameter. 

Our  model  of  the  electrolyte  solution  consists  of  a  mixture  of  charged  hard  spheres  of  radius 
q  (the  ions)  and  hard  spheres  with  a  permanent  point  dipole  of  radius  rg.  These  molecules  can  also 
have  a  sticky  dipole  [7]  interaction  which  may  preferentially  bind  the  anions  or  cations.  Analysis 
of  this  system  within  the  MSA  leads  to  explicit  expressions  for  the  thermodynamic  functions,  in 
particular,  the  Gibbs  energy  of  solvation.  Our  previous  analysis  [1-3]  was  based  on  the 
assumption  that  the  dipole-dipole  contribution  could  be  neglected,  the  Gibbs  energy  of  solvation 
being  equated  to  the  MSA  estimate  of  the  ion-dipole  contribution.  This  assumption  definitely 
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affects  the  value  of  the  polarization  parameter  estimated  from  experimental  data,  but  confirmation 
that  it  is  a  reasonable  approach  was  obtained  when  it  was  noted  [2]  that  it  leads  to  a  reasonable 
separation  of  the  components  for  cations  and  anions  when  experimental  data  for  the  alkali  metal 
halides  are  examined.  Nevenheless,  a  better  description  of  the  system  is  obtained  when  dipole- 
dipole  interactions  are  considered. 

In  the  present  manuscript,  the  model  used  previously  is  extended  to  include  consideration  of 
dipole-dipole  interactions.  The  Gibbs  energy  and  entropy  of  solvation  of  the  alkali  metal  and 
halide  ions  in  water  are  reconsidered  using  the  MSA  model  and  the  results  are  compared  with  those 
obtained  earlier. 

The  Model 

In  the  MSA  for  a  system  of  spherical  hard  ions  and  dipoles  of  different  sizes,  the  following 
expression  is  derived  [8-12]  for  the  ion-dipole  contribution  to  the  Gibbs  energy  of  solvation  in  the 
limit  of  infinite  dilution: 

gTieoTi  I,  esj  (l+^i) 

Here,  z,,  is  the  valence  of  the  ion,  e©,  the  fundamental  charge,  Es,  the  static  dielectric  constant  of 
the  solvent,  Eq,  the  permittivity  of  free  space,  rj,  the  radius  of  the  ion,  and  Nq,  the  Avogadro 
constant.  The  dimensionless  parameter  ^i,  which  depends  on  the  nature  of  the  solvent  and  the  ion, 
is  given  by 

^i=rs/(Xri)  (2) 

where  Ts  is  the  radius  of  the  solvent  molecule  represented  as  a  sphere,  and  K  the  MSA  polarization 
parameter.  In  earlier  wwk  [6-9],  X  was  calculated  from  the  dielectric  constant  of  the  pure  solvent 
Es  using  the  Wcnheim  relationship  [13] 

X2(U1)4=16Es  (3) 

Although  this  value  of  the  polarization  parameter  may  be  appropriate  for  estimating  dipole-dipole 
interaction  energy,  it  gives  poor  estimates  of  the  Gibbs  solvation  energy  when  this  is  assumed  to 


equal  Gid  [1-3]  Thus,  values  of  X  appropriate  for  cations  and  for  anions  were  determined 
separately  on  the  basis  of  experimental  data  for  the  alkali  metal  and  halide  ions,  respectively. 
Thus,  for  a  monovalent  monoatomic  cation,  Gid  is  given  by 
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Gjd(C+)  =  -  f  1  -  -"i  — 

V  (ri 


(ri+Sc) 


where 


6c  —  rg/Xc 


Xc  being  the  corresponding  cationic  polarization  parameter.  Similarly,  for  monovalent  monoatomic 


anions 


Gid(A 


^  .  No  (2ieo)2  _  K  _ 1 


EsJ  (r,+6A) 


where 


5a  -  ^sAa 


and  Xa  is  the  polarizadon  parameter  for  the  anion.  In  our  simple  theory,  6^  and  5c  are  clearly 
adjustable  parameters.  The  difference  between  them  is  due  to  the  fact  that  real  solvent  molecules 
are  non-spherical,  and  that  there  are  specific  chemical  bonds  between  the  ions  and  the  solvent 
which  are  different  for  the  cations  and  the  anions.  The  difference  could  be  interpreted  to  result 
from  an  attractive  or  repulsive  sticky  interaction,  and/or  a  different  distance  of  closest  approach,  or 


both.  Thus, 


ts  +  As 


Or,  in  other  words. 


Xi  =  X 

rs+Asi 


=  «: 


The  second  contribution  to  the  Gibbs  solution  energy  comes  from  repulsive  dipole-dipole 
interactions  [10, 11]  and  is  given  by  the  equation 


5 


(zieo)2(es-l)2 
64  jt  eoesri(l+4i)2 


£s  +  A(A±3I. 

L  2a+l)2j 


(10) 


Here,  the  parameter  X  is  assumed  to  be  that  given  by  the  Wenheim  equation,  that  is,  the  value 
appropriate  for  dipole-dipole  interactions  [10,1 1].  This  expression  may  be  simplified  considerably 
when  one  considers  the  range  of  values  typical  for  X.  For  water  whose  dielectric  constant  is  78.3 
at  25  °C,  X  is  equal  to  2.65.  The  corresponding  value  of  the  ratio  X(X+3)/2(X+l)2  which  appears 
in  the  denominator  of  the  term  in  square  brackets  in  eq.  (10)  is  0.56.  This  is  negligible  with 
respect  to  the  value  of  and  may  be  neglected.  The  ratio  appearing  in  the  numerator  of  this  term, 
namely,  (3X+2)/(X+l)  is  equal  to  2.7.  Thus,  the  expression  for  Odd  may  be  rewritten  as 


O*,  =  (I  . 


(11) 


64:teori  I,  Esj  (l+4i)^ 

The  simplified  expression  has  the  advantage  that  the  parameter  X  does  not  appear  on  its  own  but 
only  in  the  ratio  This  fact  is  important  in  applying  the  expression  to  experimental  data.  Finally, 
on  the  basis  of  previous  work  [11]  only  a  fraction  6  of  the  dipole-dipole  term  is  used  in  estimating 
the  Gibbs  energy  of  solvation  so  that  the  final  expression  for  Gs  becomes 

Gs  =  Gid  +  9  Gdd  (12) 

The  expressions  for  the  corresponding  entropies  are  obtained  from  the  temperature 
derivatives  of  the  above  Gibbs  energies.  For  the  ion-dipole  term,  the  result  for  a  cation  is 


s.(C^)  -  (1  -  i)  “1^' 


dT 


(13) 


871  eo  LEs^  ^  (ri+5c)  (ri+dc)^ 
where  d5c/dT  is  the  temperature  derivative  of  5c  which  is  obtained  from  the  temperature  derivative 
of  Xc-  Comparing  eqs.  (4)  and  (13),  it  is  easily  shown  that 


«  Gid  (C-*-)  ^  87t  eoG,d^((^)  ^ 

bid  (C  )  ^ J  =  -  ^  ^ 


(14) 


The  corresptHiding  expression  for  an  anion  is 
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Skl(A-)['l 


Girt  (A-)  ^  87t  £oGi(i^(A-)  d5A 
£j.2  dT  ■  No(zieo)2  dT 


(15) 


From  these  equations  it  follows  that  two  additional  parameters  are  required  to  estimate  the  entropy 
of  solvation  of  a  salt,  namely,  the  temperature  derivatives  of  the  polarization  parameters  for  the 
cation  and  anion. 

From  the  expression  for  Gdd,  one  obtains  the  following  equation  for  the  corresponding 
entropy: 


_2_  d^ 
dT 


Sdd  =  Gdd 

The  derivative  dX/dT  may  be  found  from  the  Wertheim  equation  from  which  one  obtains 


(16) 


^  _ 

dT  ■ 


X(l+X) 

3X+1 


1  des 


(17) 


Assuming  that  the  fraction  0  is  independent  of  temperature,  the  resulting  expression  for  the 


entropy  of  solvation  is 

Ss  =  Sid  +  0  Sdd 


(18) 


This  model  is  now  examined  with  respect  to  data  for  monoatomic  monovalent  ion  solvation  in 
water. 

Results  and  Discussion 

In  order  to  apply  the  above  model  to  experimental  data  one  must  choose  radii  for  the  ions  and 
an  extrathermodynamic  assumption  for  extracting  single  ion  solvation  energies  from  experimental 
values  for  the  salts,  namely,  the  alkali  metal  halides.  Two  sets  of  ionic  radii  give  good  fits  between 
experiment  and  theory  [2],  namely,  the  values  given  by  Pauling  [14]  and  those  extracted  from 
neutron  and  X-ray  diffracticm  experiments  [15].  The  Pauling  values  are  used  here  simply  because 
they  are  more  familiar.  The  extrathermodynamic  assumptions  used  to  separate  experimental 
quantities  into  those  for  cations  and  anions  were  examined  in  detail  for  aqueous  systems  by 
Conway  [16].  Values  of  the  Gibbs  energy  and  entropy  of  solvation  for  the  alkali  metal  cations  and 
halide  anions  estimated  by  Conway  are  summarized  in  Table  1  together  with  their  Pauling  radii. 
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In  our  previous  analysis  of  the  Gibbs  solvation  energy  [1,2]  the  contribution  of  dipole-dipole 
interactions  was  ignored  ( 0  =  0  in  eq.  ( 12)).  Analysis  of  the  data  for  the  alkali  metal  halides 
showed  that  the  best  values  of  5c  and  5a  (eqs.  (4)  and  (6))  assuming  Pauling  radii  for  the  ions 
were  82.4  and  18.0  pm,  respectively  [2].  The  quality  of  the  fit  between  theory  and  experiment  is 
really  excellent,  the  standard  deviation  for  the  calculated  values  being  1.5  kJ  mol*^  that  is, 
significantly  less  than  1  percent.  This  result  is  illustrated  in  Fig.  1  in  terms  of  a  plot  of  -Gs  against 
l/(ri-t-5i).  It  is  apparent  that  the  values  of  Gs  fit  very  well  on  a  straight  line  in  the  region  where  data 
exist  but  that  the  slope  of  the  best  straight  line  through  these  data  is  slightly  higher  (75.2  kJ  nm 
mol"')  than  that  predicted  by  theory  (69.45  kJ  nm  mol**).  As  a  result,  an  ion  with  infinite  radius  is 
predicted  to  have  a  Gibbs  solvation  energy  of  35.7  kJ  mol'*  when  this  simple  model  is  applied. 
This  can  be  attributed  to  the  failure  to  account  for  dipole-dipole  interactions  which  constitute  a 
positive  contribution  to  Gs. 

The  same  data  were  re-examined  on  the  basis  of  eq.  (12)  using  the  expression  for  Gdd 
given  by  eq.  (11).  The  fit  was  carried  out  so  that  the  resulting  estimates  of  Gid  and  Gdd  would  go 
to  zero  for  an  ion  of  infinite  radius.  The  resulting  best  values  of  5t  and  5a  were  much  smaller  than 
those  obtained  by  the  simpler  analysis,  and  are  recorded  in  Table  2.  The  corresponding  estimates 
of  Gid  are  plotted  against  l/(ri-H5i)  in  Figure  2.  As  expected  these  estimates  are  larger  in  magnitude 
than  the  experimental  values  of  Gs  also  plotted  on  the  same  graph.  Both  plots  are  linear  in  l/(r,-H5i) 
and  have  zero  intercepts  within  the  experimental  standard  deviation.  Estimates  of  0Gdd  for  the 
same  nine  ions  are  also  plotted  in  this  figure.  For  the  cations,  this  quantity  varies  finom  a  high  of 
41 .2  kJ  mol**  for  Li'*’  to  a  low  of  22.7  kJ  mol**  for  Cs'*',  and  is  approximately  eight  percent  of  the 
magnitude  of  Gid.  I®  case  of  the  anions  it  represents  a  slightly  higher  fraction  of  the  total  Gibbs 
energy  and  varies  from  a  high  of  42.7  kJ  mol**  for  F*  to  a  low  of  30.8  kJ  mol**  for  I*  ion.  The 
values  of  0Gdd  are  also  approximately  linear  in  l/(r,-^^)  with  an  intercept  corresponding  to  zero  fw 
an  infinitely  large  ion.  The  value  of  5  used  in  estimating  Gdd  *s  that  estimated  on  the  basis  of  the 
Wertheim  equation,  namely,  5 1 .6  pm.  When  the  contributions,  Gjd  and  0Gdd  are  added  one 
obtains  estimates  of  Gs  which  agree  with  those  obtained  from  experimental  data  to  within  a  few  kJ 
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mol"  I  Values  of  the  parameters  used  in  our  calculations  including  the  fraction  0  are  summarized  in 
Table  2. 


A  successful  model  for  ionic  solvation  must  also  be  able  to  estimate  correctly  the  entropy  of 
solvation.  In  order  to  make  this  calculation  one  must  have  estimates  of  the  temperature  derivatives 
of  the  parameters  5c,  5a  and  X.  The  latter  was  calculated  using  eq.(15)  which  gives  a  value  of 
-2.48  X  10"®  K"i  for  dX/dT  on  the  basis  of  the  dielectric  properties  of  water.  Assuming  rj,  the 
radius  of  a  water  molecule,  is  equal  to  137  pm  the  corresponding  value  of  d5/dT  is  0.048  pm  K"^ 
Values  of  5c  and  5a  were  determined  from  the  single  ion  estimates  of  Sg  given  by  Conway  [15]. 


In  the  case  of  cations  on  the  basis  of  eqs.  (14)  and  (18),  one  may  write 

d5c 

Y=^X 


(19) 


where 


87t  £oGid2(C^-) 

No(zieo)2 


(20) 


and 

Y  =  Ss(C+)  (^1  -  ^  +  0Sss  (21) 

Sc(C*‘)  being  the  solvation  entropy  of  the  cation.  Similar  equations  can  be  written  for  anions  using 
eqs.  (15)  and  (18).  A  plot  of  Y  against  X  for  both  ions  is  shown  in  Fig.  3.  These  data  were  fit  to 
a  straight  line  passing  through  zero  (one  parameter  fit)  as  shown,  the  resulting  slope,  which  equals 
both  d5s/dT  and  d5A/dT,  being  0.032  pm  K"^  It  is  important  to  note  that  the  present  analysis 
suggests  that  these  temperature  coefficients  are  equal  in  contrast  to  our  previous  conclusion  based 
on  an  analysis  which  ignored  dipole-dipole  interactions.  The  estimated  error  in  the  slope  is  not 
large,  being  about  4%.  It  should  also  be  noted  that  the  magnitude  of  dS^/dT  and  dSA/dT  is  about 
65  percent  of  that  found  on  the  basis  of  the  Wertheim  equation  (eq.  (17)). 

Using  the  above  coefficients  values  of  Stu  were  estimated  for  the  alkali  metal  and  halide  ions 
and  are  plotted  in  Fig.  4  together  with  estimates  of  Scu  and  experimental  values  of  8$  against  the 
experimental  Gibbs  solvation  energy  Gg.  It  is  clear  that  the  entic^y  of  solvation  is  approximately  a 
quadratic  function  of  the  Gibbs  energy  as  suggested  by  the  MSA  model. 
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Table  1.  Pauling  radii  and  standaixl  thermodynamic  parameters  of  solvation  of  monoatomic 
monovalent  ions  in  water  at  298*’  K 


Ion 

Pauling  Radius 

Gibbs  Solvation  energy 

Solvation  Entropy 

r,,  nm 

Gs/kJ  mol'i 

Ss/J  K-i  mol-i 

1. 

Li+ 

0.060 

-487.5 

-141 

2. 

Na+ 

0.095 

-387.6 

-109.7 

3. 

0.133 

-314.0 

-  74.2 

4. 

Rb+ 

0.148 

-292.7 

-  62.1 

5. 

Cs+ 

0.169 

-260.5 

-  59.1 

6. 

F* 

0.136 

-456.8 

-136.8 

7. 

a- 

0.181 

-339.7 

-  79.9 

8. 

Br 

0.195 

-325.9 

-  64.4 

9. 

I- 

0.216 

-279.5 

-  41.4 
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Legends 
Figure  1 


Figure  2 


Figure  3 


Figure  4 


for  Figures 

Plot  of  the  single  ion  Gibbs  energy  of  solvation  for  the  alkali  metal  cations  (•)  and 
halide  anions  ( ♦ )  against  the  reciprocal  of  the  ion  radius,  rj,  plus  the  MSA  parameter 
5.  The  value  of  5  for  cations  was  82.4  pm  and  for  anions,  18.0  pm.  The  solid  line 
shows  the  least  squares  fit  to  these  data,  and  the  broken  line  the  value  of  the  ion- 
dipole  contribution  to  the  Gibbs  energy  according  to  eqs.  (4)  and  (6).  The  numbers 
refer  to  individual  ions  listed  in  Table  1. 

Plot  of  the  single  ion  Gibbs  energy  of  solvation,  Gg,  the  ion-dipole  contribution,  Gid, 
and  the  dipole-dipole  contribution,  Gdd.  against  the  reciprocal  of  the  ion  radius  ri  plus 
the  MSA  parameter  6  for  the  alkali  metal  cations  (•,  o)  and  halide  anions  (4,0).  The 
value  of  5  for  cations  was  69.2  pm,  and  for  anions,  1.1  pm.  The  solid  lines  show 
the  least  squares  fit  to  each  set  of  data.  The  numbers  refer  to  individual  ions  listed  in 
Table  1. 

Plot  of  the  function,  Y,  (eq.  (17))  against  the  function  X  (eq.  (18))  for  the  alkali  metal 
cations  (•)  and  halide  anions  ( ♦ ).  The  solid  line  was  determined  by  a  one  parameter 
least  squares  fit  going  through  the  origin.  The  numbers  refer  to  individual  ions  listed 
in  Table  1. 

Plot  of  the  single  ion  entropy  of  solvation,  Sg,  the  ion-dipole  contribution,  Sid,  and 
the  dipole-dipole  contribution,  Sdd.  against  the  Gibbs  energy  of  solvation  Gg  for  the 
alkali  metal  cations  (•,  o)  and  halide  anions  (4,0).  The  curves  are  fitted  to  data  for 
Sid  and  Sdd-  The  numbers  refer  to  individual  ions  listed  in  Table  1. 
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